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Within the R-parity violating minimal supersymmetric standard model (MSSM), we use a hi- 
erarchical ansatz for the lepton-number violating trilinear Yukawa couplings by relating them to 
the corresponding Higgs- Yukawa couplings. This ansatz reduces the number of free parameters in 
the lepton-number violating sector from 36 to 6. Baryon-number violating terms are forbidden by 
imposing the discrete gauge symmetry Baryon Triality. We fit the lepton-number violating parame- 
ters to the most recent neutrino oscillation data, including the mixing angle #13 found by Daya Bay. 
We find that we obtain phenomenologically viable neutrino masses and mixings only in the case 
of normal ordered neutrino masses and that the lepton-number violating sector is unambiguously 
determined by neutrino oscillation data. We discuss the resulting collider signals for the case of a 
neutralino as well as a scalar tau lightest supersymmetric particle. We use the ATLAS searches for 
multi-jet events and large transverse missing momentum in the 0, 1 and 2 lepton channel with 7 TeV 
center-of-mass energy in order to derive exclusion limits on the parameter space of this R-parity 
violating supersymmetric model. 
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I. INTRODUCTION 

A main objective of both multi-purpose experi- 
ments ATLAS and CMS at the Large Hadron Collider 
(LHC) is the search for new physics beyond the Stan- 
dard Model (SM). Many of these extensions, in partic- 
ular supersymmetry (SUSY) [HQ, include new heavy 
colored states and a weakly interacting lightest new 
particle escaping detection. Thus the most generic 
signal among these models are several hard jets and 
large transverse missing momentum {jf> ) . ATLAS and 
CMS grouped their multi-jet and missing transverse 
momentum searches into 0, 1, 2 lepton studies 043, 
in order to be sensitive to different SUSY models and 
to avoid an overlap between these studies. Most stud- 
ies were recently updated to the full dataset of about 
5 fb _1 recorded in 2011 at a center-of-mass energy 
of 7 TeV. So far, no excess above SM expectations 
has been observed and strict bounds on any supersym- 
metric model or another relevant new physics model 
providing a similar collider signal can be derived. AT- 
LAS and CMS mainly concentrate on SUSY searches 
which are based on R-parity conserving (R p ) super- 
symmetric extensions of the SM fl5j . An equally well 
motivated scenario is a R-parity violating (/? p ) su- 
persymmetric SM [Iff, where the discrete symmetry 
baryon triality (B3) |l7( is imposed in order to avoid 
baryon-number violation and proton decay. The par- 
ticle spectrum is the same as for R p models. However, 
lepton (L-) number is violated and the lightest super- 
symmetric particle (LSP) is not stable any more. Thus 
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an alternative dark matter candidate may be needed 
such as the axino or gravitino [2ol [26| . In principle, 
any supersymmetric particle can now be the lightest 
supersymmetric particle (LSP) [27| . The LSP decays 
lead to observable effects at the LHC, which can be 
significantly different from models with R-parity con- 
servation [H,^ . Also, the L-violation causes massive 
neutrinos to emer ge i n the B3 minimal supersymmet- 
ric SM (MSSM) [19^2^ without introducing a new 
see-saw mass scale or extending the particle spectrum 
[23l |2~H . Data from neutrino experiments can be used 
to constrain the L-violating couplings (2fjj . 

Within the B3 minimal supersymmetric SM 
(MSSM), we make a hierarchical ansatz in the L- 
violating sector, relating the trilinear L-violating 
Yukawa couplings to the Higgs- Yukawa couplings, as 
first proposed in Ref. [22[. This reduces the number 
of free L-violating parameters to six. We take into 
account experimental results on neutrino oscillations, 
which amounts to five constraints (neutrino mixing an- 
gles and mass-squared differences). When additionally 
fixing the overall neutrino mass scale, this enables us 
to unambiguously determine the magnitude of the six 
L-violating parameters, removing all degrees of free- 
dom from the L-violating sector. 

Consequently, the decay properties of the LSP in 
the hierarchical B3 MSSM depend only on the exper- 
imental neutrino data. We expect no difference in the 
production and decay chains of supersymmetric parti- 
cles compared to the R p MSSM, since the magnitude 
of the L-violating couplings needs to be fairly small (of 
order 10 -5 ) in order to be in accordance with neutrino 
data. 

There have been several ATLAS and CMS searches 
as well as phenomenological studies for /? p models, 
based on resonant slepton production, multi-lepton 
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signatures or displaced vertices [3CH36l |. However, 
most of these studies constrain models where the L- 
violating couplings are either very large (for single 
slcpton production), very small (for displaced vertices) 
or where we have single coupling dominance and four 
body decays (4 lepton signature) [37| ■ Neither of these 
criteria is the case in most regions of the hierarchical 
B3 MSSM parameter space. Apart from these studies, 
the results of the ATLAS 1 lepton, multi-jet and f 
study with 1 fb _1 of data were used to restrict a bi- 
linear R-parity violating model [38| . which takes into 
account constraints from neutrino data [9|. 

In this study, we would like to re-interpret the AT- 
LAS studies with jets, f T and 0, 1 or 2 isolated leptons 

dUGl in the light of thc hierarchical B 3 MSSM. Ex- 
cept for the 2 lepton study, which uses 1 fb _1 , the stud- 
ies have been updated to 5 fb^ 1 0, QjJ, using the full 
2011 data. Since in a generic B3 MSSM, the number 
of free parameters in the SUSY breaking sector is too 
large to perform a systematic study, we work in the 
B 3 constrained MSSM (B 3 cMSSM) [HI, which im- 
poses simplifying assumptions on the scalar and gaug- 
ino masses and couplings at the unified (GUT) scale. 
It turns out that only specific regions of the cMSSM 
parameter space are phcnomcnologically viable when 
taking into account neutrino data [20j . and we focus 
on these parameter regions. As a result, there are 4 
free parameters in the SUSY breaking sector besides 
the six L- violating parameters. 

In Sect. [11], we shortly discuss how neutrino masses 
are generated in the hierarchical B3 cMSSM. We then 
describe how we fit the L-violating parameters in or- 
der to obtain the correct masses and mixing angles of 
the neutrino sector at any parameter point in the hier- 
archical B 3 cMSSM parameter space. In Sect. IIII1 we 
examine the arising collider signatures for the case of 
stau LSP and neutralino LSP scenarios. In Sect. IIV1 
we present bounds on the hierarchical B3 cSSM neu- 
trino model derived from SUSY ATLAS searches. We 
conclude in Sect. [V] 



II. HIERARCHICAL BARYON TRIALITY 
CMSSM AND MASSIVE NEUTRINOS 

A. Hierarchical Baryon Triality (B 3 ) cMSSM 

The B3 MSSM allows for additional, L-violating 
terms in the superpotential compared to the R p 
MSSM IHHH. 

W Bs = Wr, + e ab [\\ ]k L1L)E k + \' ljk L1Q)D k 

-KiLSHl]. (1) 

Li, Qi correspond to the SU(2) doublet lepton and 
quark superfields. Ei, Di are the SU(2) singlet lep- 
ton and down-type quark superfields, respectively. 
k € {1,2,3} are generation indices, a, b g {1,2} 



(ei2 = 1) are indices of the SU{2) l fundamental repre- 
sentation, while the corresponding SU(3) C indices are 
suppressed. The trilinear couplings \j k correspond to 
nine independent parameters due to the antisymmetry 
of the first two indices whereas the trilinear cou- 
plings X' i j k denote 27 independent parameters. The 
bilinear couplings 3 dimensionful couplings. 

For universal supersymmetry breaking, the bilin- 
ear L-violating couplings and the corresponding soft- 
breaking terms can be simultaneously rotated to zero 
at the unification (GUT) scale via a basis transforma- 
tion of the lepton and Higgs superfields [1^, H(| ■ How- 
ever, non- vanishing m terms (and non-aligned soft- 
breaking terms) are generated at the electroweak scale 
via the renormalization group equations [44[ ■ 

In the B 3 constrained MSSM (cMSSM), the num- 
ber of free parameters in the soft-breaking sector is 
constrained. We end up with 5 + n independent pa- 
rameters at the GUT scale (39j . 

M , M 1/2 , Aq, sgnfot), tan/?, A. (2) 

Mo, M1/2 and Aq denote the universal scalar mass, 
universal gaugino mass and universal trilinear scalar 
coupling, respectively. sgn(/z) is the sign of the super- 
potential Higgs mixing parameter and tan (3 is the ra- 
tio between the two Higgs vacuum expectation values. 
A denotes a subset of n independent dimensionless tri- 
linear L-violating couplings. 

In this work, we further restrict the number of 
free L-violating parameters: In the B3 cMSSM, the 
down-type Higgs supcrficld and the SU(2) doublet lep- 
ton superfield have the same gauge quantum numbers 
[45I I . They are indistinguishable because lepton num- 
ber is broken. Thus, the L-violating trilinear terms 
in Eq. ([T]) resemble terms in the R-parity conserving 
superpotential, 

W Rp D e ab [(y B ) ik H$L)E k + {Y D )jkH2Q h jDk\, (3) 

where (Ye)^. and (Yo)jft are the Higgs- Yukawa cou- 
plings of the lepton and thc down-type quarks, respec- 
tively. We therefore proposed the following ansatz at 
the GUT scale [22[, which can be motivated in the 
framework of Froggatt-Nielsen models (4jj 

= h ■ (Y B ) jk - li ■ {Y E ) lk , (4) 
Kjk = % ■ (Y D ) jk . (5) 

Here, £i, l\ are c- numbers. Eq. (Q| has the required 
form to maintain the anti-symmetry of the \j k in 
the first two indices. Assuming a specific form of the 
Higgs Yukawa couplings, the number of L-violating 
parameters reduces to six complex numbers. We have 
given our ansatz in the weak-current basis. How- 
ever, after EW symmetry breaking, we must rotate to 
the mass-eigenstate basis. Experimentally, only the 
PMNS and the CKM matrix are known Hi, S3- The 
explicit lepton and quark mixing matrices are there- 
fore not fully determined. In thc following, we assume 
that the lepton Higgs- Yukawa matrix is diagonal Thus, 
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we assume mixing only in the neutrino sector for the 
leptonic sector. In the quark sector, we assume left- 
right symmetric mixing. Additionally, we work in the 
limit where the down-type Higgs- Yukawa matrix is 
diagonal whereas the up-type is non-diagonal. Hence 
our specific form of the Higgs- Yukawa couplings im- 
plies mixing only in the up-type-sector. In Ref. [2l| , it 
was shown that the choice of quark mixing (e.g. mix- 
ing in the up-type versus mixing in the down-type- 
sector) does not significantly influence the numerical 
results at the low energy scale. 



B. B3 neutrino masses 

Since lepton number is violated, the neutrinos mix 
with the neutralinos, resulting in a 7x7 neutralino- 
neutrino mass matrix of rank 5. As a result, we obtain 
one massive neutrino at tree level 1391 . 



167ra G uT ^ i=1 ( Ui 



v d - 



M 1/2 



Here v d , v u and v% denote the vacuum expectation 
values of the H d , H u and sneutrino fields. However, 
experimental neutrino oscillation data suggests that 
we need at least two massive neutrinos. Since there is 
only one massive neutrino at tree-level, higher-order 
corrections need to be taken into account. Full 1— 
loop corrections to the neutrino-neutralino mass ma- 
trix have been discussed in Ref. [2l| . A good estimate 
of the size of these radiative corrections is given by 
the slepton-lepton and down-type quark-squark loop 
contribution, which are proportional to [48j 
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(7) 
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The proportionality of the loop contributions to the 
exchanged SM fermion mass in the loop further in- 
creases the effect that trilinear couplings with indices 
i33 are dominant over all other indices ijk, as is clear 
from the hierarchical ansatz in Eqs. ((4]) and ([5]). 

Ref. [ITJl noted that in large regions of cMSSM pa- 
rameter space the ratio between the tree-level neu- 
trino mass and the radiative contributions is too large 
too yield a phenomenologically viable neutrino mass 
hierarchy. However, due to RGE effects in the run- 
ning of L- violating parameters, the tree-level neutrino 
mass has a global minimum at 
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2Af 1/2 , 
M 1/2 



(9) 
(10) 



for non-zero A^ fe or Aijk, respectively. We choose Aq 
such that it minimizes the A' contribution to neutrino 
masses [Eq. Q], as explained in more detail in the 



next paragraph. Thus, in the hierarchical B3 cMSSM 
a set of 10 free parameters, 



M1/2, M , sgn(n), tan/3, £j, 4 



(11) 



fixes the full B 3 cMSSM. 

As described in Ref. [2l| , it is possible to obtain the 
experimentally measured neutrino mass squared dif- 
ferences and mixing angles by independently generat- 
ing each neutrino mass with a set of three L-violating 
free parameters. This means that 6 or 9 independent 
couplings are necessary in order to obtain the full spec- 
trum with either two or three massive neutrinos. How- 
ever, in the case of neutrinos in normal hierarchy mass 
ordering with a massless lightest neutrino, it turns out 
that one can do with only 2 couplings to explain the 
heaviest neutrino mass, m„ 3 , cf. Ref. [2l|. This is 
fortunate, because due to our hierarchical ansatz only 
£[, £\ and £2 have a significant impact on the neutrino 
sector whereas £3 generates only a negligible contribu- 
tion to the neutrino masses if it is of the same order 
of magnitude as the other couplings [7l[. Therefore, 
(6) we generate m„ 3 at tree-level via the A^k couplings, 
which are in turn determined by £\ and £ 2 . The second 
neutrino mass, m„ 2 is generated via A' i j k (determined 
by the £■) at one-loop level, whereas the lightest neu- 
trino must remain massless, m Ul « 0. 

In summary, we have 5 free L- violating parameters 
which control the neutrino sector, £' i and £1, £2- These 

* j v-i i re- 



can be used to generate non-zero m U2 and r, 
spectively, in accordance with the two mass squared 
difference and three mixing angles from experiment. 
It is not easily possible to obtain inverse hierarchy 
or degenerate neutrino masses in the hierarchical B3 
cMSSM unless £3 becomes several orders of magnitude 
larger than the other L-violating parameters. 



C. Experimental neutrino oscillation data 

Assuming three active oscillating neutrinos, the best 
global fit values of the neutrino masses and mixing 
parameters at la C.L. are given by 0, [HFJ], 

sin 2 [6»i 2 ] = 0.31 ±0.02, 
sin 2 [6» 23 ] = 0.51 ±0.06, 
sin 2 [26»i 3 ] = 0.09 ±0.02, 

Am^ = 7.59 ± 0.2 x 10" 5 cV 2 , 



Si 



-2.34 ±0.1 x 10~ 3 eV 2 
2.45 ±0.1 x 10" 3 eV 2 



where 



Am, 2 , = m 2 



(12) 



(13) 



m Ui denote the neutrino masses in order of largest 
electron-neutrino admixture. There are two large mix- 
ing angles 6*i 2 and 623- Deviating from Ref [4Sj . we use 
in Eq. (TT2"]) for #13 the best fit value recently measured 
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by Daya Bay and RENO [5fJ, GLy]. The neutrino os- 
cillation data implies at least two non-vanishing neu- 
trino masses m Vi . In this work, we usually consider 
the so-called Normal Hierarchy (NH) scenario, where 
0. 



Amj^ > and m 




D. Numerical results 

For each cMSSM point we fit the L-violating pa- 
rameters £i and £\ to the best-fit Normal Hierarchy 
neutrino mass data in Eq. (|12[) . Wc perform this fit 
by minimizing the x 2 function 



(14) 



where f° s are the central values of the N b s experi- 
mental observables in Eq. (fT2j). y^ oftsus y are the cor- 
responding numerical predictions and Si are the la 
uncertainties. We calculate the low energy mass spec- 
trum and couplings with S0FTSUSY3 . 3 . lH] . The nu- 
merical minimization of our \ 2 function is done with 
the program package MINUIT2 [53j]. Details of our nu- 
merical procedure can be found in Ref. [2l1 | . Here, we 
present an example solution where we translate the 
best fit values £i and £\ into the corresponding values 
of the trilinear L-violating couplings at the unification 
scale: 



A 133 


= 1.72 
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A233 


= 2.74 
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= 1.13 
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A 233 


= 3.89 
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A 333 


= 3.11 


10 


-5 



(15) 

Wc have used M = 100 GeV, M 1/2 = 500 GeV, 

tan/3 = 25, sgn(^) and Aq X ' « 2A/ 1 / 2 - As one can 
see, the A^33 and A -33 couplings are between 0(1O~ 5 ) 
and 0(1O~ 6 ). All remaining trilinear L-violating cou- 
plings are at least one order of magnitude smaller, 
below O(10 -7 ). The couplings A'233 and A 333 tend 
to be the largest trilinear L-violating couplings. In 
Fig. [TJ we display the best fit value of A' 23 3 in the 
M -M 1 / 2 plane. We see that the magnitude of the L- 
violating couplings does not strongly depend on M 
and Mx 1 1- Furthermore, the relative magnitude of the 
L-violating couplings to each other remains roughly 
the same throughout the parameter space. 

Recall that the parameter £3 is not fixed by the neu- 
trino oscillation data in the normal hierarchy scenario. 
However, we assume that £3 is of the same order of 
magnitude as l\ and £2, setting £3 = £2 in the rest of 
our paper [72j |. 

We have checked all low energy constraints on the L— 
violating trilinear couplings [54, 55 1. However, in our 
case the couplings are too small to have an observable 
impact on any low energy observables. 
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FIG. 1. Best-fit values of the L-violating coupling A233 at 
the unification scale in the Mo~M 1 / 2 plane, fixing Ag X ' ~ 
2M 1/2 , tan/3 = 25 and sgn(p) = +1. 



III. COLLIDER SIGNATURES 

In this section, we investigate possible collider sig- 
natures of the hierarchical B3 cMSSM at the LHC. 
The best-fit values of the L-violating couplings to neu- 
trino data are too small to have an observable effect 
on the resonant production of supersymmetric parti- 
cles. Thus, pair production of colored sparticles via 
strong interactions is the dominant production channel 
at the LHC. Only if slcptons and gauginos are much 
lighter than the colored sparticles, their production 
rate becomes comparable. The produced sparticles 
cascade decay into the LSP. In our parameter space, 
we can have either a stau LSP or a neutralino LSP [73| . 
The final state collider signature is determined by the 
decay properties of the LSP candidate. In the B3 
cMSSM, the LSP is almost always short-lived and de- 
cays within the detector via the L-violating interac- 
tions [74|]. We now describe the final state signatures 
of stau LSP and neutralino LSP scenarios separately 
after describing the numerical tools used. Then we go 
on to discuss in which regions of Mq — M\j2 parameter 
space they occur. 



A. Numerical tools 

The low energy mass spectrum and couplings 
are calculated with S0FTSUSY3 . 3 (52|. The decay 
widths of the relevant sparticles are obtained with 
IsaJet7.64 [H] and IsaWigl.200. However, the de- 
cay channels of the neutralino LSP via the sneutrino 
vevs and the Ki term are absent in IsaWigl.200. 
Therefore, we calculate decays via the bilinear L— 
violating couplings with SPrieno3 . 1 [57} . We combined 
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all decay widths in order to calculate the branching ra- 
tio of the sparticles. We use the parton distribution 
functions MRST2007 LO modified [Hi. Our signal 
events are generated with Herwig6 . 510 59J. The cross 
sections are normalized with the NLO calculations 
from Prospino2.1 [60l | assuming equal renormaliza- 
tion and factorization scale. Our events are stored in 
the Monte Carlo event record format StdHep5.6.1. 
We take into account detector effects by using the fast 
detector simulation Delphesl . 9 [6l|, where we choose 
the default ATLAS-likc detector settings. Our event 
samples are then analyzed with the program package 
ROOT [H| and we calculate the 95% and 68% confi- 
dence levels (CL) of the exclusion limits with TRolkc 



B. Stau LSP decay 

In the parameter region where the lighter stau fi 
is the LSP, pair produced squarks and gluinos at the 
LHC cascade decay into the LSP, producing jets and 
taus (tau-neutrinos) along the way, 



pp qq/qg/gg nn + 2j + x, 



(16) 



where j and X denote jets and additional particles 
of the process (such as r or v T ) 1 respectively. Note 
that we can have more than 2 jets in the final state 
if the process involves gluinos. These additional jets 
are included in X , which wc discuss in more detail in 
Sect IIIIDI For example, right-handed squarks decay 
into a jet and the lightest neutralino, which then typ- 
ically decays into a stau and a tau with a branching 
ratio of one, 



qn -> jjXiXi 



J J TT nn 



(17) 



The stau then directly decays into two SM fermions 
via the trilinear L-violating couplings A133, A233 and 
A 3jfc , cf. Fig [3] Decays via the Ai33 couplings are 

is 



dominant, even though the decay width via A 3 i . 
enhanced by a factor of Nq = 3 and the A 3ji . cou- 
plings are generally larger. However, the lightest stau 
is mostly right-handed and thus the coupling of the 
stau via A' is suppressed due to the small admixture 
with the left-handed stau. Additionally, the stau de- 
cay via A333 into a top and bottom quark is kinemat- 
ically forbidden or suppressed in large regions of pa- 
rameter space. Stau decays via A 311 and A 322 are heav- 
ily suppressed due to the smallness of the couplings. 

In principle, the stau can also mix with the charged 
Higgs boson via K3 and decay via the two-body decay 
mode f — > tv. However, we have numerically checked 
that stau decays via bilinear operators are always sub- 
dominant in our model. We define a 

• benchmark point BP1 in the stau LSP region 
with 

A/ = 100 GcV, A/ 1/2 = 500 GcV, tan ,3 = 25, 




FIG. 2. Schematic characterization of the stau LSP decay 
in the hierarchical B3 cMSSM. 



This benchmark point is characterized by lightest neu- 
tralino, lighter stau, gluino and squark masses of 205 
GeV, 162 GcV, 1146 GcV and 1012 GeV, respectively. 
The dominant LSP branching ratios for BP1 are given 
by 



Br(f 1 ^ t v c ) 


= 0.26 


Br(ff -> T~Vy) 


= 0.21 


Br(f 1 _ — > c~v T ) 


= 0.26 


Br(f-f — > jit~i/ T ) 


= 0.21 


Br(f 1 ~ — > s c ) 


= 0.04 



(18) 

Note that the branching ratios into different decay 
channels are roughly independent of the stau mass as 
long as the final state masses are negligible. 

Roughly half of the staus decay into a charged lep- 
ton and neutrino, the other half decays into a tau and 
neutrino. Note that we only denote electrons or muons 
as leptons in this paper. Since one third of taus decays 
leptonically, we expect final state collider signatures 
with either 0, 1 or 2 leptons from the decaying stau 
LSPs, for 12%, 46% and 42% of events, respectively: 

0£ + 2v + 2T had + 2j + X 
U + 2(4)u+lT had + 2j + X 

2£ + 2(4, 6)1/ + 2j + X (19) 

where £ denotes an electron or muon and Thad denotes 
a hadronically decaying tau. If the lcpton[s] in the 1£ 
or 21 channel come from a leptonically decaying tau, 
the number of neutrinos increases from 2 to 4 [6], as 
shown in brackets in Eq. (|19p . Due to the Majorana 
character of the neutralino, both neutralinos can decay 
into like-charged staus and hence we can have same- 
sign leptons in the final state. 



C. Neutralino LSP decay 

In the hierarchical B3 cMSSM, the lightest neu- 
tralino eigenstate is generally bino-like. The produc- 
tion process is given by 



PP 



qq/qg/gg -► Xixl + 2j + x . 



(20) 



sgn(//) = 1 and A, 



(A') 



2M 1 



/2 



The neutralino LSP can either decay via a trilinear 
L-violating operator into three SM fermions or via 
neutralino-neutrino mixing (proportional to the bilin- 
ear L-violating couplings and the sneutrino vevs) into 
a gaugc/Higgs boson and a lepton, cf. Fig. [31 



G 




\ b W(Z) 

FIG. 3. Schematic characterization of the three-body (left) 
and two-body (right) decay modes of the neutralino LSP 
in the hierarchical B 3 cMSSM. 



For relatively low sfermion masses in the propaga- 
tor, the trilinear three-body decay modes dominate 
because the bilinear L-violating couplings are only 
generated radiativcly via RGE running and the sneu- 
trino vevs are determined to be relatively small from 
radiative electroweak symmetry breaking. However, 
in parameter regions with heavy sfermions, the bilin- 
ear two-body decay mode becomes dominant because 
the three body decay mode suffers from phase space 
suppression and heavy virtual sfermions in the propa- 
gator. 

First, we discuss the case where the lightest neu- 
tralino dominantly decays via the trilinear LNV cou- 
plings, for which we define 

• benchmark point BP2 with 

M = 200 GeV, M 1/2 = 400 GeV, tan/3 = 25, 
sgn(/i) = 1 and w 2M 1/2 

This benchmark point is characterized by lightest neu- 
tralino, lighter stau, gaugino and squark masses of 
163 GeV, 213 GeV, 937 GeV and 846 GeV, respec- 
tively. We obtain the following LSP branching ratios 
for BP2: 

Br(x? -> %bb) = 0.31 

Br(x° -> %bb) = 0.20 
Br(x° -> W ± ^ =F ) = 0.21 
Br(x? -> W*^) = 0.05 

Br(x? -> = 0.13 

Br(x? -»■ trh°) = 0.08 (21) 

The branching ratio of the three-body decay modes 
(the \i ~~ * vbb channel) is roughly 51%. However, 
for this benchmark point the two-body L-violating 
decays via bilinear L-violating couplings already have 
a sizable contribution to the LSP decays. The electron 
(electron-neutrino) channel is suppressed compared to 
the muon decay channel because A' 133 ~ 0.3A 233 , cf. 
Eq. (fl~5j) . Therefore, about 90% of our leptons are 
muons. Summing up the various decay channels and 
including the gauge boson branching ratios, roughly 
72% of neutralinos decay without leptons, 19% with 



one lepton and 7% with two leptons. This leads to 
52%, 27% and 14% of events with 0, 1 and 2 leptons 
from LSP decays, respectively. 

Assuming the cascade decay processes of Eq. (|20p. 
dominant final state signatures are then given by 

0£ + 2p + 2bb + 2j + X 
U + lv+ bb + W had + 2j + X 

2i + 2v + bb + 2j + X (22) 

Next, we discuss the decay properties of the light- 
est neutralino in a region where the two-body decays 
dominate, 

• benchmark point BP3 with 

M = 600 GeV, M 1/2 = 400 GeV, tan/3 = 25, 

sgn(/i) = 1 and A V) sa 2M 1/2 

The lightest neutralino, lighter stau, gluino and squark 
masses of BP2 are 164 GeV, 579 GeV, 961 GeV and 
1010 GeV, respectively. Here, the LSP decay chan- 
nels are the same as for BP2; however, the branching 
ratios differ drastically: 



Br(x? 


-> vibb) 


= 0.04 


Br(x? 


-> v T bb) 


= 0.03 




= 0.40 


Br(x? - 


> W ± T T ) 


= 0.14 


Br(x? 




= 0.27 


Br(x? 


-> %h°) 


= 0.12 



(23) 



Since here the scalar masses (Mo) are fairly large, 
the two-body neutralino decay modes via bilinear 
L-violating couplings or sneutrino vevs dominate, 
amounting to 93%. Therefore, there are only half as 
many neutralinos decaying into the 0£ channel as for 
BP2; twice as many decay into the \l and 21 chan- 
nel. This results in final state signatures with 0,1 or 2 
leptons at 24, 37 and 27%, respectively. Typical final 
state signatures are given by 

0£ + 2v + 2Z° ad/vv + 2j + X 

U + lv+ Z° ad/vv + W had + 2j + X 

2l + 2v+ Z° ad/vv +2j + X (24) 

As mentioned before, the electron decay channel is 
suppressed by roughly a factor of 10 compared to the 
muon decay channel. Additionally to the channels 
mentioned in Eq. (f2"4")l. there are 12% of events with 3 
or 4 leptons from LSP decay. 

D. Scan in the Afo — M1/2 plane and kinematical 
distributions 

In the subsequent numerical analysis, we perform 
a scan in the Ma — M1/2 plane. For this, we define 
a benchmark region (BR) which contains the three 
benchmark points defined above (BP1, BP2, BP3): 
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FIG. 4. The iso curves show the logarithmic ratio between 
three-body and two-body decay modes of the neutralino 
LSP in our benchmark region. In the stau LSP region, the 
two-body stau decay modes via the trilinear RPV cou- 
plings are always dominant. 

• Benchmark region BR (where Mo, Mi/ 2 
free): 

tan/? = 25, sgn(/i) = 1 and A^' ] w 2M 1/2 

BP1, BP2 and BP3 each lie in distinct sections of 
the BR: stau LSP region, neutralino LSP region domi- 
nated by three-body decays and neutralino LSP region 
dominated by two-body decays, respectively. This is 
depicted in Fig 01 where the ratio between three- and 
two-body decay modes of the neutralino LSP is dis- 
played. The two-body Xi decay modes dominate at 
large M 1 / 2 and Mq. As one can also see in this figure, 
the stau LSP region within our BR is approximately 
given by 

M 1/2 > 3 M - 80 GeV, (25) 

since the lightest neutralino mass is driven to larger 
values by the large Mi/ 2 . In general, the lighter stau 
mass eigenstate is mostly right-handed. 



In Sect. Ill Bl we discussed that the absolute magnitude 
of the L-violating parameters as well as the relative 
magnitude between them does not vary significantly 
with Mo and M]/ 2 . This implies that the LSP de- 
cay branching ratios are hardly affected by variations 
of the L-violating parameters within our BR. How- 
ever, the decay modes are importantly affected by two 
points, as illustrated in Fig. [4] 

(A) Whether we are in the stau or neutralino LSP 
region 



(B) The ratio between three- and two-body decay 
modes within the neutralino LSP region. 

In the stau LSP region, the 1 and 2 lepton channels 
are dominant for large regions of parameter space. The 
lepton channel only becomes significant once the stau 
becomes heavier than the top-quark. Then, hadronic 
stau decays via A^ 33 contribute significantly and the I 
and 2 lepton studies perform much worse, resulting in 
a "cutoff" of the sensitive region for stau masses above 
the top mass. Now, the lepton channel could further 
exclude parameter space; however, since this region 
extends well above M 1 / 2 ~ 500 GeV, we expect that 
the amount of data collected is not yet large enough 
to make exclusion possible. In the neutralino LSP re- 
gion dominated by three-body decays, we expect the 
lepton channel to be the best, whereas in the case of 
two-body decays, the 2 lepton channel should perform 
better. 

We now come to a discussion of possible additions 
to the final state particles from "AT" [as contained in 
Eqs. (HHJ) and (|2"D)) ] and the most important distribu- 
tions for our benchmark region. 

Additional jets can arise from gluinos in the hard 
process, since the gluino decays into quark and (vir- 
tual) squark, leading to more jets in the final state [75j. 
Besides gluino pair and gluino-squark production, 
gluinos can occur in squark decays if Mi/ 2 <C Mo. 
For example, in BP3 the gluinos are lighter than the 
squarks and a sizable fraction of the squarks decay 
into a gluino and a quark which then decays via vir- 
tual squark and quark. Thus, we expect a higher jet 
multiplicity than for BP1 or BP2, where rriq < rrig. 
This is illustrated in Fig. [5] (i). There, we show the 
distribution of the number of jets for our three bench- 
mark points as well as for a i? p -conserving version 
of BP2 and BP3 with a stable LSP (denoted "BP2 
RPC" and "BP3 RPC", respectively). One can see 
that for BP2 RPC, there are on average only 2-3 
jets because here squarks typically decay into a neu- 
tralino/chargino and a quark, whereas for BP3 RPC, 
there are 3-4 jets. Comparing BP2 RPC to BP2, we 
expect up to 4 additional b-jets from the neutralino 
LSP decays [Eq. (|22p], and thus the distribution peaks 
around Aj ot =5 — 6, cf. Fig. [5] (i). Similar observa- 
tions can be made for BP3. Here, there are more jets 
from the (R-parity conserving) decay chain involving 
gluinos. However, on average there are less jets from 
neutralino LSP decays, Eq. (f2~4")l . such that the distri- 
bution also peaks at Aj Ct = 5 — 6. In the stau LSP 
case (BP1), the distribution peaks at Aj Ct = 3 — 4. 
Here there arc only few jets which can be attributed 
to X (ie. gluino decays), as discussed above. 

Further leptons in the final state can emerge in the 
cascade decays of the SU(2) doublet squarks. The lat- 
ter decay into charginos and neutralinos with domi- 
nant SU(2) gaugino composition, which are typically 
xf and x.2 m the cMSSM. xt and x 2 subsequently 
decay either into slepton and lepton or gauge bo- 
son/Higgs and the lightest neutralino. However, this 



leads to isolated leptons in only ~ 15% of events in our 
case, as is illustrated in Fig. [5] (ii) by the Ng distribu- 
tions for BP2 RPC and BP3 RPC. The reason for 
this is that in BP2, the f\ is much lighter than the 
other sleptons, whereas the latter are heavier than \2 
and xt- Thus X2 an( i xt dominantly decay into fr 
and fv, respectively. About one third of these r's de- 
cay leptonically, leading to final state leptons. In BP3, 
all sleptons are heavier than xt and X2 an( i hence the 
latter preferably decay into a gauge/Higgs boson and 
the lightest neutralino. Comparing BP2 RPC and 
BP3 RPC with the corresponding /J scenarios, we 
clearly see that there are significantly more leptons for 
BP2 and BP3 due to lcptonic decays of Xi- However, 
there are more entries in the lepton bin for BP2 and 
BP3 than expected from Eqs. (|22[1 and (f2~4"]l . because 
some of the leptons are non-isolated or too soft or do 
not fall into the acceptance region of the tracking sys- 
tem. The same holds for BP1, which has overall the 
largest number of isolated leptons; nevertheless the ra- 
tio between events with 1 lepton and leptons is still 
less than predicted from Eq. (|T9|) . 

In Fig.[5](iii), we present the missing transverse mo- 
mentum distribution. Here, we clearly see that BP1 
has the hardest distribution among all R p violating 
distributions. Note that for the two other R p violating 
scenarios the missing transverse energy distribution is 
much softer compared to the respective R p conserving 
scenarios, due to the LSP decays. 



IV. NUMERICAL RESULTS: EXCLUSION 
LIMITS ON HIERARCHICAL B 3 CMSSM 
PARAMETER SPACE 

In this section, we further constrain the hierarchical 
B3 cMSSM parameter space using data from the LHC 
at yfs = 7 TeV with an integrated luminosity of up to 
5 fb _1 . We focus on recent ATLAS studies with 0,1 or 
2 isolated leptons, several jets and large missing trans- 
verse momentum. A short overview over the ATLAS 
studies used is given in Table [U Full details of objects 
reconstruction, definitions of all kinematical observ- 
ablcs and event selection cuts of all three analyses can 
be found in the resp ective ATLAS publications 0- 
H (0 lepton), @, [lof (1 lepton) and [3 (2 leptons). 
We have chosen these analyses because they only rely 
on simple objects such as electrons, muons, jets and 
missing transverse momentum in the final state. Thus, 
we do not rely on complicated tau reconstruction and 
b-tagging algorithms, which are difficult to simulate 
with the detector simulation Delphesl .9 [6l[. In par- 
ticular, difficulties arise in reconstructing hadronically 
decaying taus [28|. Also, the published ATLAS studies 
for supersymmetry involving taus [64[ or b-jets [(35| in 
the final states have smaller cross-sections or smaller 
efficiencies than the multi-jet, large p and lepton 
searches. Thus, we expect the "simple" 0-2 lepton 
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FIG. 5. We depict (i) the number of jets iVjet, (ii) the num- 
ber of isolated leptons Ni with pr > 20 GeV and (iii) the 
missing transverse momentum ("ETMISS") for our bench- 
mark points BP1, BP2 and BP3. Additionally we dis- 
play an R p version of BP2 and BP3 ("BP2 RPC","BP3 
RPC"), where the neutralino LSP is kept stable. We gen- 
erated 40000 events for each benchmark point. 



analyses to perform better with the current amount of 
data. So far, the experimental data is in agreement 
with the SM background expectations. We use their 
results in order to derive the 68% and 95% CL ex- 
clusion regions in the Mq-M^ii parameter space. We 
plan to investigate exclusion limits arising from third 
generation studies and multi-lepton studies in a future 
publication. 





Olept-SRE-m 


llept-3j 


21ept-OS-4j 


N t 





1 


2 




6 


3 


> 4 


T 
Pjots 


>(130, 60, 60, 60, 40, 40) >(100, 25,25) >(100, 70, 70, 70) 


t T 


> 160 


> 250 


> 100 




> 1200 


> 1200 




t T 

m cff 


> 0.15 


> 0.3 




C 


4.7 fb _1 


4.7 fb" 1 


1.0 fb _1 



TABLE I. The main cuts used in the ATLAS studies used 
in this collider study. More details concerning the cuts can 
be found in the relevant ATLAS studies (0 lepton [jj, 1 
lepton [l^l and 2 lepton f]2])- Nt denotes the number of 
isolated leptons, Nj e t the number of jets and pf et s specifies 
the minimal transverse momentum which is required for 
these jets. ft T gives the minimal value of missing transverse 

momentum of the event, m^ c ' the minimal (inclusive) ef- 
fective mass and £ denotes the total integrated luminosity 
at 7 TeV. 



ATLAS and CMS have recently published confer- 
ence notes which found that the lightest Higgs is at 
least heavier than 117.5 GeV at 95% CL [66|, |67|. In 
the hierarchical B3 cMSSM, the lightest Higgs is typ- 
ically rather lighter than 116 GeV, because the value 
of Aq is necessarily fixed to be positive and similar in 
magnitude to 2M 1 / 2 , cf. Sect. Ill Bl This means that 
the stop mixing cannot become very large and thus 
the loop contributions to the lightest Higgs mass are 
moderate. We have checked various values of tan/3 
and both sgn(/i) = ±1; however, we found that the 
Higgs mass does not become larger than 117 GeV for 
Mo, M 1 / 2 < 1 TeV. Therefore, the exclusion limits 
derived from this lightest Higgs mass bound would by 
far exceed the exclusion limits derived from the 0, 1 
and 2 lepton channels mentioned above. However, it 
could be possible to soften the bound if we extend the 
field content of the hierarchical B3 MSSM by a singlet, 
i.e. working in the next-to minimal SSM (NMSSM) 
[68l - i70| . We leave this topic for a future investigation 
at a time when there is more certainty regarding the 
lightest Higgs mass. 

Before applying the model independent cross sec- 
tion limits from the ATLAS searches to our neutrino 
model, we checked that the Monte Carlo tools are cor- 
rectly tuned. Therefore, we generated 20000 events 
for each grid point in the Mq-Mi/ 2 plane in the R- 
parity conserving cMSSM. We determined the 95% CL 
exclusion region in the Mq-M 1 / 2 plane for the AT- 
LAS "Hepton-3j" study (cf. Table Q} and verified that 
our results are compatible with the interpretation from 
ATLAS within ±30 GeV. We now discuss the 0, 1 and 
2 lepton channels in detail. 



A. lepton channel 

ATLAS has used the lepton channel as one of the 
first search channels for supersymmetry 0-13 • So far, 
they have collected a total luminosity of about 4.7 fb _1 
at the center of mass energy of ^fs = 7 TeV. From the 
non-observation of an excess, we can derive exclusion 
limits on the hierarchical B 3 cMSSM. The ATLAS 
lepton channel is divided into several signal regions 
(SR). For all signal regions, the cut on p and the 
minimum requirement on pj et of the first two most- 
energetic jets are identical. However, the number of 
jets and the minimum pj ct cut for the remaining jets 
as well as the cut on m*g c and on the ratio p T /m e s 
differ for the different signal regions. 

We have examined all signal regions after apply- 
ing the object reconstruction described in their study 
and found that we obtain the strictest exclusion lim- 
its for the "SRE-m" signal region, which demands 6 
jets, m 1 ^ 1 > 1200 GeV and ^ > 0.15, cf. Table H 
We show the resulting plot in the Mq-Mi/2 plane 
in Fig. [5] The exclusion limit peaks at Mo ~ 200 
GeV. This is the region where the neutralino LSP de- 
cays dominantly via three-body decays Xi ~^ v bb, 
c. f. Fig. 2J It was to be expected that the "SRE- 
m" signal region gives good exclusion limits for this 
type of scenario, because if both neutralinos decay via 
Xi — > vbb, we expect at least 6 parton level jets (in- 
cluding b-jets). Also, we have only moderate be- 
cause of the three-body decay of the neutralino, and 
therefore more events survive in the "SRE-m" than 
in the "SRE-t" scenario (where m*| c > 1500 GeV). 
Finally, leptons from the cascade decays of SU(2) dou- 
blet squarks into Xi an d Xo are suppressed, since the 
latter dominantly decay into xt ~~ ^ ™ ano - X2 ~ * 

For increasing Mo, the exclusion region decreases 
to lower Mx/2 values. We can see in Fig. 2] that the 
two-body decay mode of the neutralino becomes more 
important here. Thus, an increasing number of the 
neutralino LSPs decay into a gauge boson and a lep- 
ton and less &-jets are expected in the final state, so 
that less events pass the kinematical cuts on the final 
state jets. Another effect is that for larger Mo, the 
production cross section decreases. 

Directly to the left of the peak at M « 200 GeV, 
the limit drops off sharply because here the LSP be- 
comes the fi and there are significantly less events 
with 6 jets and no leptons. However, Mx/ 2 < 350 
GeV can still be excluded at 95% CL. We would like 
to point out that in principle, it is possible to obtain 
better exclusion limits (up to Mi/ 2 < 400 GeV) in the 
stau LSP case by using a signal region with only 4 
or 5 jets. However, the 1 lepton study performs even 
better and therefore we go not into detail about the 
results from these signal regions here. 

We do not consider the region with a LSP lifetime 
exceeding cr = 15 mm, since the ATLAS searches for 
supersymmetry require prompt LSP decays. In Fig. 
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FIG. 6. Exclusion limit on our benchmark region, where 

f a' ) 

tan/3 = 25, sgn(/i) = 1 and Ag ~ 2A{ 1 / 2 , from the iso- 
lated leptons, 6-jets and MET ( "Olept-SRE-m" ) ATLAS 
study. The white region is excluded at 95% confidence level 
(CL), the light blue is excluded at 68% CL. The grey lines 
denote the squark masses, the dashed black lines denote 
the gluino masses (each in GeV). The black line delineates 
the region (below) where the lifetime of the LSP becomes 
larger than cr > 15 mm. In this region, the exclusion 
limits are not applicable because the ATLAS study rejects 
leptons and jets which do not originate from the primary 
vertex. 

EJ the region below the solid black curve highlights a 
LSP with a lifetime ct > 15 mm. 



of the llept-3j signal region are listed in Table Q] and 
the resulting plot is shown in Fig. [7] Almost half of the 
events in the stau LSP region decay into final states 
with 1 lepton, cf. Sect. IIIIBI Note also that the 1 lep- 
ton study [lOj demands the most stringent cut on p 
among the 0, 1 and 2 lepton studies. In the stau LSP 
region with direct (two-body) leptonic decays, much 
more missing transverse momentum is produced than 
in the neutralino LSP region. In particular in the neu- 
tralino LSP region with dominant three-body decays 
into ubb, the amount of jj> T is greatly reduced compared 
to the stau LSP region. Moreover, much less charged 
leptons arise from the neutralino decay Additional 
leptons from the cascade decays are also heavily sup- 
pressed. Thus, we have a sharp drop of the acceptance 
in the crossover region between the stau and neutralino 
LSP region. For larger Mo values, eventually the two- 
body neutralino decay modes become dominant over 
the three-body decay mode. However, the hard cut 
on p T still rejects many signal events in this region. 

Note that the ATLAS signal region with 1 lepton 
and 4-jets is also sensitive to the neutralino LSP re- 
gion besides the stau LSP region. This explains why 
in the old 4-jet signal region with 1 fb _1 in the muon 
channel, ATLAS was able to constrain the bilinear ft p 
model presented in Ref. Q (with two-body neutralino 
decays) quite well. However, having in mind that in 
our case we have addtional three-body decays and in 
the new 5 fb _1 study, the cuts are more stringent cuts 
than the 1 fb _1 version and not optimized for our type 
of scenario, the resulting exclusion limits on the neu- 
tralino LSP region are weaker than the limits derived 
in the 2 lepton channel as shown below. 



B. 1 lepton channel 

Refs. H [l(| search for multi-jet events with large 
missing transverse momentum and exactly one iso- 
lated lepton. Similarly to the lepton channel in the 
previous subsection, the 1 lepton channel was one of 
the first supersymmetry search channels and the cur- 
rent integrated luminosity is 4.7 fb _1 at the center of 
mass energy of 7 TeV. They consider signal regions 
with 3- or 4-jets with different kinematic configura- 
tions, which are optimized for the R p cMSSM with a 
large mass difference between the gluino and the LSP. 
Additionally, they include a soft-lepton signal region 
which is sensitive to scenarios with small mass split- 
ting between the sparticlcs. 

Comparing the results for the different signal re- 
gions, we observe that the 3-jet signal region ("llept- 
3j" ) provides us with the best overall exclusion limits 
in the stau LSP region up to Mm ~ 500 GeV (i.e. 
better than the limits from any other signal region in 
the to 2 lepton channels). The main kinematic cuts 



C. 2 lepton channel 

The ATLAS study based on final states with two 
leptons and missing transverse momentum (l2| has not 
yet been updated to include more than 1 fb _1 of data. 
The search is divided into opposite-sign (OS), same- 
sign (SS) and flavour-subtraction (FS) signal regions 
where up to 4 jets are demanded besides exactly 2 lep- 
tons and a cut on p T . We find that we obtain the best 
exclusion limits with the OS signal regions. The three 
OS regions differ in the ft T cut, the number of jets and 
the corresponding minimal pj cts cut. As in the case 
of the 1 lepton channel, the OS studies with the hard- 
est transverse missing momentum cut ("21ept-OS-2j" , 
f T > 250 GeV) are quite sensitive to the stau LSP re- 
gion where two staus decay leptonically. However, in 
the 2 lepton channel the obtained exclusion limits are 
~ 50 GeV weaker than in the "llept-3j" study. This 
is due to the stringent cuts on m"^ c and on the ratio 
p T /m e ff in the "Hcpt-3j" search channel, which yield 
better signal isolation and background suppression. 

The OS and 4-jet channel with a moderate f T cut of 
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FIG. 7. Exclusion limit on our benchmark region, where 
tan/3 = 25, sgn(/x) = 1 and Ag X ' « 2Mi/ 2 , from the 
1 isolated lepton, 3-jets and MET ("Uept-3j") ATLAS 
study The white region is excluded at 95% CL, the 

light blue is excluded at 68% CL. The grey lines denote the 
squark masses, the dashed black lines denote the gluino 
masses (each in GeV) . The black line delineates the region 
(below) where the lifetime of the LSP becomes larger than 
CT ^ 15 mm. In this region, the exclusion limits are not 
applicable because the ATLAS study rejects leptons and 
jets which do not originate from the primary vertex. 




Mo(GeV) 



FIG. 8. Exclusion limit on our benchmark region, where 
tan/3 = 25, sgn(/x) = 1 and Aq ' m 2Mx/ 2 , from the 2 iso- 
lated opposite-sign leptons, 4-jets and MET ("21ept-OS- 
4j") ATLAS study. The white region is excluded at 95% 
CL, the light blue is excluded at 68% CL. The grey lines 
denote the squark masses, the dashed black lines denote 
the gluino masses (each in GeV) . The black line delineates 
the region (below) where the lifetime of the LSP becomes 
larger than cr > 15 mm. In this region, the exclusion 
limits are not applicable because the ATLAS study rejects 
leptons and jets which do not originate from the primary 
vertex. 



100 GEV ("21ept-OS-4j"), described in Table U pro- 
vides us with the best exclusion limits for Mo > 300 
GeV, where the neutralino LSP decays dominantly via 
two-body decays as shown in Fig. [SI We notice a slight 
dip for smaller Mo (Mo ~ 200 GeV), where there are 
dominant three-body neutralino decays. Here, as dis- 
cussed in the previous subsections, parton-level lep- 
tons from the neutralino LSP decays or from the cas- 
cade decays of the SU(2) doublet squarks are heavily 
suppressed and the exclusion limits from the lepton 
channel are more stringent. For even smaller values of 
Mo, we are in the stau LSP region and the exclusion 
limits improve again. However, as discussed in the last 
paragraph, the cuts are not optimized for a stau LSP 
scenario. The E™ lss cut is the weakest among all three 
analyses in Table U and the kinematic requirements on 
the jets are harder compared to the "llept-3j" search 
channel. 

For Mo 3> Mi/2, the gluino is generally lighter than 
the squarks and thus we expect a higher jet multiplic- 
ity and in general more jets passing the kinematic cuts. 
However, much less transverse momentum is generated 
compared to the R-parity conserving case or the stau 
LSP region. Thus, the "21ept-OS-4j" yields the better 
overall exclusion region in the neutralino LSP region 
with dominant bilinear RPV decays due to the softer 
E™ iss cut compared to " Olept-SREm" . One further 
remark on the number of leptons in the final state: for 



Mi/2 ^ -^0, the SU(2) doublet squarks decay via a 
wino-like gaugino is quite sizable, although we have 
the competing decay channel via an off-shell gluino. 
These wino-likc gauginos again dominantly decay into 
gauge bosons providing additional leptons in the final 
state. 



V. SUMMARY AND CONCLUSION 

We introduced a hierarchical ansatz for the L- 
violating trilincar Yukawa couplings in the B3 cMSSM. 
Here, the trilinear LNV Yukawa couplings are related 
to the Higgs Yukawa couplings via six independent pa- 
rameters £i and t\. We have then determined the best 
fit values of the £i and l\ in order to obtain phenomeno- 
logically viable neutrino masses and mixing angles. 
It is possible to quasi unambiguously determine the 
L-violating sector as well as the value of the SUSY 
breaking scalar coupling Aq from neutrino oscillation 
data. We discussed the final collider signatures in the 
stau LSP and neutralino LSP scenarios at the LHC 
and finally used the ATLAS searches in jets and large 
missing transverse momentum with 0, 1 and 2 isolated 
leptons in order to find the 95% and 68% CL exclu- 
sion limits in the Mq-M 1 / 2 plane for fixed sgn(/x) and 
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tan/3. We can exclude squark masses below 800 GeV, 
and gluino masses below 700 GeV (for squark masses 
below 1 TeV) at 95%. These limits become more strin- 
gent at 68% CL, by roughly 100 GeV. Compared to 
the case of the R-parity conserving cMSSM, we obtain 
weaker limits because generally we have more jets and 
leptons and less p T due to the LSP decays. 

We want to conclude with a short discussion of how 
we can improve a future collider study for our model 
or similiar R-parity violating models. There are a 
number of studies in which R-parity violating collider 
signatures are investigated, as mentioned in the in- 
troduction. Many of these studies consider multilcp- 
ton (Ni > 3) signatures in association with much less 
missing transverse energy than in our study. They typ- 
ically assume, however, a single non-zero \jk without 
third generation indices, ie. k £ {1, 2}, so that the 
number of lepton is enhanced. In our model, the LSP 
dominantly decays via Aj33 or \' i33 couplings involving 
third generation decay products, or via neutralino- 
neutrino mixing involving gauge boson decay prod- 
ucts. However, the branching ratio of the LSP into 
leptons is still considerably large (between 19% and 
47%) and therefore the lepton multiplicity is higher 
than in R-parity conserving models. Also, the aver- 
age pj distribution of the signal leptons will be rel- 
atively hard due to the large phase space of the two 
body decay channels of the LSP into SM fcrmions. 
For example, in BP2 the hardest lepton has on aver- 
age pj = 80 GeV. In BP2 RPC, the hardest lepton 
has a mean value of pj = 60 GeV. Demanding one 
(two) lepton(s) with moderate pj cuts might be ad- 
vantageous to isolate the signal. As an alternative, we 
can also apply a kinematical cut on the scalar sum of 
all the leptons' pj . 

Decays via trilinear couplings with third generation 
indices are dominant in large regions of parameter 
space in our model. Therefore, we expect a substantial 
proportion of events with third generation SM parti- 
cles in this parameter region. For example, we ex- 
pect a large number of taus and b-jets in BC1 and 
BC2, respectively. In BP2, about 50% of all events 
have at least one b-jet. This is in sharp contrast to 
BP2 RPC where only 13% of all events have a b- 
jet in the final state. Requiring hadronically decaying 
taus or b-jcts in the final state should help to sup- 
press the SM background. However, for the parameter 
region around BP3, the LSP dominantly decays via 



neutralino-neutrino mixing. Here, we do not expect 
third generation particles in the final state in abun- 
dance. 

The increase in jet and lepton multiplicities due to 
LSP decays in our model happens at the cost of less 
missing transverse momentum compared to the R— 
parity conserving case. For example, in BP3 RPC 
we have on average p = 213 GeV because the stable 
neutralino LSP escapes detection. In BP3 we obtain 
a mean value of p = 123 GeV due to neutrinos from 
the LSP decay. In many studies the effective mass, 

Meg =f T +J2pf ets , (26) 

is used to "measure" the effective SUSY mass scale. 
However, they assume a stable LSP and thus M e g re- 
ceives a sizable contribution from p . Our signatures 
tend to look softer than those of most R-parity con- 
serving scenarios because some of the decay products 
of the LSP are not included in the sum in Eq. (|26p. 
A useful discriminating variable to increase the signif- 
icance of our signal could be the scalar sum of missing 
transverse momentum, all jets, leptons and hadronic 
taus, 

s t = t T + E pfets +E^+E pL. • ( 27 ) 

For example, the ratio of Eq. and Eq. (J2ZJ) is 0.85 
for BP3. 

Finally, it is difficult to constrain the region Mi / 2 < 
230 GeV in our model due to the finite lifetime of the 
LSP, since many supcrsymmetry searches only recon- 
struct leptons and jets which originate from the pri- 
mary vertex. We thus conclude that allowing events 
with displaced vertices would certainly be advanta- 
geous to establish bounds in the low M\/2 region. 
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